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Abstract: ln order to analyze skull variation in the genus Macaca , seventy-seven cranicdental variables were taken 
from eleyen species. They were first defined seven functional units comprising three anatomical regions . Twenty-seven vari- 
ables were finally selected to carry out the morphology of the whole skull. The data, organized in these ways, were examined 
to discover variations between and within the various species. The methods used were Principal Components Analysis 
(PCA) and Discriminant Function Analysis (DFA) . PCAs of the functional units anatomical regions, and the whole skull 
provided similar, though not identical, separations of species clusters in both sexes separately . These differences in structure 
could be related to size. sexual dimorphism, diet , ecology, classification and phylogeny. The question of size should have 
been easy to settle . Unfortunately, this is not the case. In this study where the raw data are measurements of the specimens, 
the main differences should be size. However „the size differences seem to occur in both the first and second (independent) 
multivariate axes. ln some analyses the size differences between the species are biggest and appear in the first axis ,ln other 
analyses it is the separation between the sexes (and these too are largely size ) that are the biggest and appear in the first 
axis. Yet in other analyses, both of these size separations, though still orthogonal to one another „present in the combination 
of the first two axes. This certainly implies that a single axis of body size is not present and that shape differences have not 
been isolated form size differences. lt also implies that sexual dimorphism is a complex matter. As a result ,the question of 
the relationships between the species is therefore also complex. One cluster of species that includes M . fascicularis , 
M. sinica and M . radiata was significantly isolated from all others regardless of level of analysis. This relationship is quite 
different from that proposed on the anatomy of the reproductive organs (Delson , 1980; Fooden, 1976 . 1980). 
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There have not been many craniodental morpho- 
metric studies of macaques { Pan, 1998; Pan et al., 
1998). Such results as there are suggested that the 
studies provide an excellent anatomical system with 
which to test hypotheses of morphological integration 
and classification. Some results were alse useful in re- 
vealing the relationship between craniodental morpholo- 
gy and its underlying developmental bases” . Questions 
still remain to be tested. To what extent can phenotypic 
results provide evidence clarifying controversies about 
phylogeny? To what degree can variation between 
species be separated in association with functional 


adaptation and phylogenetic inertia? 
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Phylogenetically, macaques have been categorized 
into four species-groups but with different definitions 
based on the morphology and structure of the reproduc- 
tive organs (Delson, 1980; Fooden, 1976, 1980) . These 
concepts have been tested by results from other re- 
search fields, for instance, fossil records and external 
features (Eudey, 1980), biological and karyological 
studies (Fooden et al., 1989; Melnick et al., 1985; 
Zhang et al., 1990), and behaviour (Chan, 1996; 
Fooden, 1975). 

Regarding differences in functional adaptation, a 
detailed comparison of diets and ecology of macaque 
species was provided by DeClue et al. (1992). The 
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various species have been described as adapted to a 
wide variety of habitats and living conditions, and as 
displaying an extensive range of diets and feeding be- 
havior. These variations might be relevant to corre- 
sponding species differences in morphologic structure of 
the teeth and cranium. A limited study of this type has 
been undertaken by Hylander (1979) and Takahashi et 
al .( 1994). 

The study aim is to provide further tests using 
craniodental morphometric analysis with a large number 
of variables and species. This study attempts to do this 
by examining craniodental variation at three distinct 
levels of anatomical organization, individual functional 
units, anatomical regions and the whole skull . Implica- 
tions for functional adaptation and phylogenetic rela- 


tionship will be examined . 
1 Materials and Methods 


Eleven species of the genus Macaca (Table 1) 
were used in this study. The materials are housed in a 
number of institutes, universities and museums in vari- 
ous parts of the world (see Acknowledgments ). All 
specimens were from adult macaques as judged by the 
full eruption of M3s.They were of wild shot animals or 
of animals that had died of natural causes. 

Seventy-seven linear dental, mandibular and cra- 
nial variables (Table 2) were measured on each speci- 
men using digital calipers. Length measurements were 
made parallel to the median sagittal plane, and breadth 


{or width ) measurements were made in the coronal 


Table 1 Species and number (by sex) used in this study 








Species Male Female Unknown Total 
M . mulatta 20 20 40 
M fuscata 11 11 22 
M . sinica 14 8 22 
M . assamensis 20 20 40 
M . radiata 12 12 I4 
M . arctoides 20 20 40 
AM. sylvanus 11 10 21 
M . nemestrina 20 20 40 
M . negra 8 8 16 
M . fascicularis 20 20 40 
M . thibetana 11 8 l 20 


Total 167 157 1 325 





(Pan and Oxnard, in preparation 


plane. Heights were measured at the right angles to 
length and breadth. 

Principal Components Analysis (PCA) was ap- 
plied to the resulting data. Clusters of major contribut- 
ing variables were identified using, as a significance 
criterion, minimum loadings of 0.300 (Hair, et al.. 
1992}. The variables were studied in three successive 
steps to focus on different anatomical aspects. In the 
first step variables were grouped into 7 units on the ba- 
sis of function . Some variables were selected and divid- 
ed into three anatomic regions secondly . Twenty-seven 
variables were finally used to analyze the morphology of 
the whole skull. This division of skulls into three inves- 
tigational levels; functional units, anatomical regions 
and the skull as a whole, was aimed at investigating the 
ways in which the different parts with different fanc- 
tions separate the sexes and the species” ( Pan, 1998). 

Since a big sample size from l1 species of the 
genus Macaca was used in this study .it will be difficult 
in figuring out the separation between species if all in- 
dividuals were plotted on diagrams of PCA and DFA, 
So only means of loading scores of each sex or species 


were shown. 
2 Results 


Eigenvalues and eigenvectors for each analysis of 
each individual sex were provided, Those for the sexes 
taken together are not listed since both the name and 
number of the variables making significant contribu- 
tions to the first two axes were very similar either to fe- 
males or males taken separately . 

2.1 Principal Components Analysis 
2.1.1 Functional units 

Seven units, the mandible, lower teeth, upper 
teeth, lower face ,upper face, cranium and the calvaria, 
were analyzed separately by PCA. In these initial com- 
parisons, the patterns of separations of the species and 
the sexes are used heuristically in order to increase un- 
derstanding of the effects of analyzing the different 
anatomical regions . 


@Females The separations between species de- 
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Table 2 Percentage of total variation, eigenvalues and eigenvectors for the first two PCA axes for variables of functional 
units in the female macaques 











Mandible Lower teeth Upper teeth 
i 2 1 2 1 2 
Eigenvalue 10.94 1.46 Eigenvalue 9.64 1.66 Eigenvalue 10.68 1.33 
Percentage 63.4 9.1 Percentage 60.3 10.4 Percentage 66.7 3.3 
Cum. per 63.4 TES Cum. per 60.3 70.6 Cum, Per 66.7 75.0 
Figenvectors Figenvectors Eigenvectors 
LBCB 0.682 0.619 Li1MDL 0.249 0.921 LIIMDL 0.458 0.637 
LBMB 0.936 0.089 LI2MDL 0.562 0.707 LI2MDL 0.647 0.423 
LIAW 0.430 0.791 LCMDL 0.509 -0.310 CCMDL 0.769 0.147 
CONDYLL 0.549 -0.2331 LF3MDL 0.701 -0.211 UP3MDL 0.851 0.040 
CONDYLW 0.384 ~ 0,201 LP4MDL 0.388 0.099 UP4MDL 0.885 -0.156 
CONMI 0.956 ~0, 143 LMIMDL 0.875 0.166 LMIMDL 0.859 - 0.214 
CONM3 0.877 -0.221 LM2MDL 0.929 -0.034 LM2MDL 0.901 - 0.194 
MAMI 0.863 ~ 0,296 LMIMDL 0.909 - 0.186 UM3MDL 0.891 - 0.196 
MAT 0.725 ~ 0. 132 L11 BLL 0.623 0 O02 LI1BLL 0.620 0.569 
MANDLAL 0.896 0.011 LIZBLL 0.741 - 0.071 UIŻBLL 0.764 0.310 
MANDH 0.899 -0.070 LCBLL 0.544 0.160 UCBLL 0.814 - 0.020 
MANDSYM 0.896 0.172 LP3BLL 0.845 0.066 UP3BLL 0.867 0.043 
BICONDY 0.955 -0, 144 LP4BLL 0.914 - 0.091 UP4BLL 0.966 - 9.056 
MANCORYW I 0.863 -0.068 LM1BLLM 0.900 - 0.062 UMIBLLM 0.395 -0.263 
SYMWID 0.384 0,055 LM2BLLM 0.935 - 0.136 UM2BLLM 0.926 -0.214 
LIRLALY 0.728 0,296 LM3BLLM 0.896 -0.148 UM3BLLM 0.904 -0.136 
Maxilla Cranium Calvaria Upper face 
I 1 I 2 I 2 I 2 

Eigenvalue 4.64 0.72 Eigenvalue 6.10 0.96 Eigenvalue 4.30 0.98 Eigenvalue 3.64 0,71 
Percentage 77.4 II-9 Percentage 76.3 12.1 Percentage 61.4 14.0 Percentage 70.6 8.9 
Cum. Per 77.4 39.3 Com. Per 76.3 33.3 Cum. Per 61.4 75.5 Cum, Per 70.6 99.5 
Eigenvectors Eigenyectors Eigenvectors Eigenvectors 
MUZL 0.926 -0.I07 BCRANL 0.951 -0.028 CALVL 0.926 -0.117 MUZL 0.863 - 0.039 
PALLENG 0.959 -0.109 CRANL 0.980 0.0IS BPORW 0.911 -0.117 BIORBW 0.837 0.181 
PALWID 0.905 -0.143 CRANW 0.99 0162 MIDPARW 0.895 -0.040 INTORBW 0.726 0.633 
UBCB 0.914 -0.070 BIZYGW 0.934 -0.062 POSTORB 0.618 0.397 PIRH 0.861 -0.027 
UBMB 0.933 -0.081 POSTORB 0.557 0.68%  OCCH 0.594 -0.663 PIRW 0.755 -0.476 
UAW 0.582 0.812 OCCH 0.567 -0.674 FORMAGL 0.550 0.590 UFACEH 0.848 - 0.197 

ANTBASI 0.976 -0.041 FORMAGW 0.880 0.085 INFRMAL 0.505 0.062 

BIZYGW 0.908 - 0.092 


The variables reaching significant contribution level are in bold. For the definition of variables please see Pan, 1998; Pan, et al. , 1998. 


fine primarily three clusters and are mainly seen along 
the first axis (Fig.1). Except for some overlap between 
the second and the third clusters in the upper face, this 
finding is clear. (The number of species in the studies 
involving units in the mandible is one less than in the 
other regions because of missing data for M . nigra ).To 
facilitate comparisons, three clusters of species are out- 
lined by convex polygons in the diagrams. These poly- 
gons are not assumed ,at this point,to have any biologi- 
cal meaning with regard to the species clustering, but 
rather as providing a guide to the differential patterns 
produced by the different anatomical areas. In other 
words ,the polygons are here a heuristic device for re- 
vealing similarities in the results in the different levels 
of analysis. The first cluster ( polygon) includes the 
species M . fascicularis , M . radiata and M . sinica ; the 


second, the species M.assamensis, M . nemestrina , 
M. mulatta and M nigra; and the third. the species 
M , arctoides , M . thibetana , M . fuscata and M . sylvan- 
us. 

The first two axes in each study account for more 
than 70.0% of the total variation (Table 2) .In PC1 on- 
ly one variable, LIIMDL, does not show a significant 
positive contribution. In the second axis the variables 
showing significant positive contributions are LBCB, Ll- 
AW, LIIMDL, LI2MDL, LI1MDL, UI2MDL, UIIBLL, 
UIZBLL, UIAW, POSTORB, FORMAGL and IN- 
TORBW . Those with significant negative contributions 
are LCMDL,OCCH and PIRW. 

Males Again. in general terms the same three 
clusters of species are evident .In contrast with the situ- 


ation in females, there is somewhat greater overlap a- 
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mong these clusters in four units; lower teeth, lower 
face, cranium and the calvaria. In general, the separa- 
tion between clusters is similar to that in females al- 
though individual separations are smaller (Fig.2}. 

In males, more than 66.7% of the total variation is 
accounted by the first two axes. As with females, all 
variables in PC1 show significant positive contributions 
(Table 3). The variables showing significant positive 
contributions to PC2 include LBCB, LIAW,LTRLALY, 
LIIMDL, LIZMDL, LCMDL, LP3MDL, LIIBLL, 
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Plots of the first two principal components of the analyses of individual functional units in female macaques 


LCBLL, UI1MDL, ULAW, OCCH, FORMAGL. FOR- 
MAGW, UCMDL, UCBLL, BIORBW and INTORBW. 
Those with significant negative contributions are 
MAM1, LM2MDL, LM3MDL, LM1BLLM . LM2BLLM, 
UBMB, POSTORB, OCCH, MUZL, PIRW and U- 
FACEH. 

Two sexes together The analysis of the sexes 
taken together are illustrated in Fig.3.Three patterns of 
variation were found. The first pattern is in the 


mandible, lower teeth and the upper teeth, in which the 
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Fig.2 Plots of the first two principal components for functional units iu male macaques 
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Table 3 Percentage of total variation. eigenvalues and eigenvectors for the first two PCA axes for variables of functional 
units in the male macaques 




















Mandible Lower teeth Upper teeth 
1 2 I 2 1 2 
Figenvalue 10.79 1.27 Eigenvalue 9.31 1.94 Eigenvalue 9.66 1.81 
Percentage 67.5 7.9 Percentage $8.2 12.1 Percentage 60.4 11.3 
Cum. per 67 5 75.4 Cum. per 58.2 7.3 Cum. Per 60 4 71.1 
Eigenvectors Egenvectors Eigen ectors 
LBCB 0.779 0.505 LIIMDL 0.439 0.660 LIIMDL 0.484 0.602 
LBMB 0.914 0.048 LEMDL 0.624 0.394 URMDL 0.736 -0.013 
LIAW 0.599 0.700 LCMDL 0.601 0.304 UCMDL 0.590 0.534 
CONDYLL 0.618 - 0.040 LP3MDL 0.598 0.531 UP3MDL 0.645 0.046 
CONDYLW 0.825 -0.147 LP4MDL 0.886 0.019 UP4MDL 0.874 -0.105 
CONM1 0.940 -0.163 LM1MDL 0.829 0.244 UM1MDL 0.867 -0.184 
CONM3 0.865 -0.252 LM2MDL 0.887 - 0.303 UM2MDL 0.905 -0.207 
MAMI 0.745 - 0.335 LM3MDL 0.888 - 0.305 UM3MDL 0.399 - 0.180 
MAT 0.715 -0.240 LIIBLL 0.509 0.473 UITBLL 0.484 0.691 
MANDLAL 0.961 0.002 LI2BLL 0.772 D 258 UI2BLL 0.779 0.241 
MANDH 0.897 -0.071 LCBLL 0.562 0,334 UCBLL 0.597 0.458 
MANDSYM 0.865 0.273 LP3BLL 0.864 0.049 UP3BLL 0.854 0.056 
BICONDY 0.922 - 0.182 LP4BLL 0.868 -0.179 UP4BLL 0.361 - 0.058 
MANCORW1 0.800 -0.163 LM1BLLM 0.850 - 0.343 UMIBLLM 0.876 - 0.310 
SYMWID 0.848 - 9.028 LM2BLLM 0.890 - 0.313 UM2BLLM 0.884 -0.317 
LTRLALY 0.739 0.311 LM3BLLM 0.881 -0.282 UM3BLLM 0.874 -0.34 
Maxilla Cranium Calvaria Upper face 
I 2 I 3 1 2 I 2 
Eigenvalue 4.74 0.52 Eigenvalne 6.28 0.73 Eigenvalue 3.68 1.00 Eigeavalue 5.54 0.92 
Percentage 78.9 8.7 Percentage 78.5 9.1 Percentage 52.5 [4,2 Percentage 69.2 1.5 
Cum. Per 78.9 87.6 Cum. Per 78.5 87.6 Cum Per $2.5 66.7 Cnm. Per 69.2 80.8 
Engenvectors Eigeavectors Eigenvectors Eigenvectors 
MUZL 0.919 0.081 BCRANL 0.953 -0.097 CALVL 0.826 -0.227 MUZL 0.770 -0.400 
PALLENG 0.932 -0.023 CRANL 0.9309 0.104 BPORW 0.913 -0.176 BIORBW 0.826 0.332 
PALWID 0.891 -0.278 CRANW 0.949 -0.039 MIDPARW 0.849 0.094 INTORBW 0.701 0.619 
UBCB 0.917 0.205 BIZYCW 0.926 0.022 POSTORB 0.762 -0.086 PIRH 0.876 0.016 
UBMB 0.852 -0.429 POSTORB 0.714 -0.407 OCCH 0.518 -0.489 PIRW 0.761 -0.313 
UIAW 0.812 0.459 OCCH 0.617 0.733 FORMAGL 0.601 0.419 UFACEH 0.851 -0.366 
ANTBASI 0.964 0 023 FORMAG™W 0.482 0.694 INFRMAL 0.922 0.186 
BIZYEW 0,923 -0.038 





The variables reaching significant contribution level are in bold. 


separation between species is mainly revealed in PCI 
and that within species (i.e. between sexes) in PC2. 
The second pattern, almost opposite to the first one, is in 
the lower face and the upper face. The third is in the 
cranium and the calvaria, in which the variation. both 
within and between species, is oblique, though still or- 
thogonal in both PC] and PC2. 
2.1.2 Anatomical regions 

The aggregation of the functional units into 
anatomical regions comprised the masticatory apparatus 
(mandible, lower teeth and the upper teeth), the face 
(lower face and upper face) and the cranium as a whole 
(cranium and calvaria).In each of these analyses, all 
variables in PCI make significant contributions (Table 
4). 


Two notes are important here. The original inten- 


tion was to examine data based upon the upper jaw 
(maxilla) with the expectation that this would present a 
picture like that in the other sections of the masticatory 
apparatus. When, however. the raw data for the upper 
jaw (maxilla) were examined more critically, it was re- 
alised that they were not measures of the upper jaw but 
of the lower face . It is for this reason that the upper jaw 
or maxilla does not appear as a separate region. Again, 
though it might be expected that the pattern of descrip- 
tion might go from masticatory apparatus through face to 
cranium,the face was found to be intermediate between 
the masticatory apparatus and the cranium (as might 
have been expected). It is therefore easier to describe 
the two extremes first and to place the face last. 
‘DMasticatory apparatus Twenty-six variables se- 
lected from the mandible , lower teeth and the upper 
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Table 4 Percentage of the total variation eigenvalues and eigenvectors for the first two principal components in 


the three anatomic regions 





Masucatory apparatus 











Female Male 

1 2 1 2 
Eigenvalue 15.50 2.53 15.02 57.80 
Percentage 59.6 9.7 2.71 10.4 
Cum. Per 59.6 69.4 57.8 68.2 
Eigenvectors 
LBCB 0.736 0.059 0.792 0.258 
LBMB 0.914 -0.227 0.888 -0.219 
LIAW 0.615 0.453 0.666 0.504 
CONM1 0.868 -0.379 0.853 -0.164 
MAMI 0.742 - 0.468 0.689 -0.211 
MANDLALY 0.887 - 0.330 0.883 -0.02 
MANDH 0.819 -0.312 0.881 -0.015 
BICONDYL 0.864 -0.355 0.883 -0.21t 
ULIMDL 0.531 0.606 0.606 0.608 
UI2MDL 0.719 0.278 0.744 - 0.126 
UCMDL 0.793 0.729 














Female Male 
1 2 1 2 

UP4M BL 

UM1MDL 0.831 - 0.028 0.862 -0.275 
CN BLL 0.633 0.258 0.487 0.630 
CEBLL 0.773 0.153 0.760 0.067 
CCBLL 0.767 - 0.103 0.656 0 249 
LP4BLL 0.872 - 0.012 0.831 - 0.154 
UM1BLLM 0.898 - 0.094 0.823 - 0.394 
LILMDL 0.348 0.303 0.476 0.734 
LIZMDL 0.630 0.505 0.626 0.266 
LCMDL 0,525 -0.226 0.597 Q. 168 
LP4MDL 0.922 0.065 0.861 - 0.090 
LM1MDL 0.856 0.130 0.812 -0.234 
LCBLL 0.543 0.204 0.560 0.263 
LP4BLL 0.386 0.041 0,861 -0.202 


LM1BLLM 


Face as a whole 














Cranium as a whole 
Female Male 

1 2 1 2 
Eigenvalue 7.79 1.20 6.56 1.38 
Percentage 59.9 9.2 50.5 10.60 
Cum. per 59.9 69.1 50.5 61.1 
Eigenvectors 
CALYL 0.953 - 0.053 0.872 0.175 
BPORW 0.925 - 0.104 0.910 - 0.099 
MIDPARW 0.865 0.006 0.781 - 0.386 
INTFRMAL 0.887 - 0.243 0.920 0.097 
POSTORB 0.595 0.598 0.667 -0.197 
TEMFOSL 0.864 - 0,058 0.338 0.006 
TEMFOSW 0.842 -0.038 0.762 -0.002 
ZYGH 0.703 - 0.038 0.650 0.305 
OCCH 0.579 — 0.431 0, 530 0.335 
FORMAGL 0.496 0.049 0.498 - 0.246 
FORMAGW 0.924 —0. 180 0.367 - 0.623 
GLENOL D.495 0.722 0.303 0.685 
GLNOW 0.834 0.152 0.804 0.059 


The vanalles reaching significant contribution level are in bald. 


teeth were used to carry out this analysis based on the 
criteria described above. 

a.Females The first two axes account for 69.4% 
of the iotal variation.The same three clusters of species 
as in the analyses of separate functional units were 
found (Fig.4).The first axis mainly disperses the clus- 
ters in order of species size (the species with the small- 
esi and the largest body size are allocaied on the most 
left and right parts of PC1,separately).The firs) cluster 
is discrete but ihe second one partially overlaps with the 
species in the third cluster. The second axis also ar- 
ranges the species in order of size. Overall, an interme- 
diaie , diagonal, axis would best separate the three clus- 


ters in the first two axes. 





Female Male 

1 2 1 2 
Eigenvalue 7.30 1.00 7.07 1.33 
Fercentage 66.4 9.1 64.3 12.1 
Cum. per 66.4 75.5 64.3 76.3 
Eigenvectors 
MUZL 0.525 0.050 0.893 -0.237 
BIORBW 0.808 -0.270 0.720 0.528 
INTORBW 0.676 -0.340 0.562 0.639 
PALLENG 0.994 0.106 0.919 -0.212 
PALWID 0.864 0. 166 0.861 - 0.205 
UBCB 0.870 0.204 0.852 - 0.351 
UBMB 0.928 0.091 0.900 0.093 
UAW 0.493 0.727 0.699 - 0.478 
PIRH 0.832 - 0.303 0.814 0.345 
PIRW 0.731 -0.253 0.698 0.128 
UFACEH 0.830 0.030 0.822 0.028 


In the second axis, CONM1,MAM1,MANDLALY. 
MANDH and BICONDYL, make significant negative 
contributions, in contrast to four variables, LIAW, 
UI1MDL, LI1MDL and LI2MDL, which have significant 
positive eigenvectors (Table 4) .In other words, the an- 
terior teeth play a particular role in the separations in 
the second axis. 

b.Males The first two components accouni for 
68.2% of the total variance (Table 4). The separation 
among the three clusters in the firsi axis is not as clear 
as in the females (Fig.4) although, again, the firsi clus- 
ier is well separated from the other two. However, quite 
differently from females, the separation between the 


clusters is ambiguous in the second axis as there is al- 
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Fig.3 Differences between sexes and species on the first two principal components of the analyses of functional units 
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Fewer variables than in the females make signifi- LIJMDL,UIJMDL and LI2MDL. All are anterior dental 
cant contributions to the second axis. As with females. dimensions as in females. The variables making nega- 


however,the variable showing the greatest eigenvectors tive contributions have relatively low eigenvectors, only 
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two of the dental dimensions, UMIBLLM and 
LM1BLLM, are significant. 
c.Females and males 


species in PC1 is almost solely due fo the body size of 


The separation between 


species (Fig.5)}.The difference between the sexes with- 
in species is virtually totally in the second axis. These 
patterns are very similar to those obtained from the vari- 
ous functional units. 

@Cranium as a whole Thirteen variables were 
chosen from the cranium and the calvaria for this analy- 
sis based on the criteria already presented . 

a.Females The eigenvalues and eigenvectors 
forthese analyses are listed in Table 4,The first two ax- 
es account for 69.1% of the total variation. The first 
cluster is yet again well separated from the other two 
which overlap partly. The dispersion of species in the 
second axis is similar to that of the calvaria analysis. 
That is, the three clusters overlap completely . 

In PC2 the positive eigenvectors ,the breadth of the 


postorbital constriction (POSTORB) and glenoid length 





0.0 
PCI 


Cranium as a whole 





05 15 
PCI 


WAAT. ERRANA SRA (EX) 


O00 0D http:/www.cqvip.com 


315 


(GLENOL), are contrasted mainly with the negative 
eigenvector, the occipital height (OCCH). 

b. Males 
61.19% of the total variation (Table 4). The arrange- 


ment of the clusters is somewhat similar to that of the 


The first two components constitute 


cranium (Fig.2 and 4); the three species in the first 
cluster are significantly separated from the rest of 
species and there is a strong overlap between the sec- 
ond and third clusters . 

The variables making significant positive contribu- 
tions to the second component are ZYGH,OCCH and 
glenoid length (GLENOL) and foramen magnum 
breadth (FORMAGW ). INFRMAL and FORMAGW 
make significant negative contributions to this axis. 

c.Females and males The separation between 
sexes is in the first axis, and that between species is in 
the second axis (Fig.5). The three species, 
M.radiata, M.sinica and M./fascicularis are very 
closely grouped and significantly separated from the rest 


of the species along the second axis. Macaca thibetana 
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shows a unique position in the second axis. 

(Face as a whole Eleven variables were select- 
ed from the lower and the upper face for this analysis. 

a.Females The eigenvalues and eigenvectors are 
listed in Table 4. The first two components account for 
75.5% of the total variance. The separation and the ar- 
rangement of the clusters are rather different from those 
of the lower face, but somewhat similar to those of the 
upper face (Fig. 1 and 4).The three clusters all overlap 
totally in PC1. Some separation between clusters is 
found in the second axis where the first cluster is sepa- 
rated from the other two which strongly overlap 
(Fig.4). 

The variables showing significant negative contri- 
butions to the second axis are the interorbital breadth 
{INTORBW ) and the height of the piriform aperture 
(PIRH). Only one variable has a significant positive 
contribution, the maxillary incisor alveolar breadth 
(UIAW). 

b.Males The first iwo components account for 
76.3% of the total variation. Overall, the arrangement 
of the clusters in the two axes is rather different from 
that of either of the individual units (Fig.2 and 4) .The 
dispersions of the clusters in the first axis , however, are 
very similar to that of the lower face; there is strong 
overlap between all three clusters .In PC2 the first clus- 


ter is separated from the other two which strongly over- 
lap. 

Biorbital breadth ( BIORBW ), interorbital breadth 
(INTORBW ) and PIRH have significant positive eigen- 
vectors and are contrasted with UBCB and alveolar 
breadth (LIAW) which make significant negative con- 
tributions to the second axis. These variables, as in fe- 
males ,are mainly linked with the anterior breadths. 
The difference by sex with- 


in species lies mainly in the first component (Fig.5), 


c.Females and males 


that between species, however, is in the second axis. 
Both of these differences are primarily related to size 
(the main difference between both species and sexes ) 
but in this case, the effect of size between sexes is the 
larger. This pattern is quite similar to those of the indi- 
vidual units (Fig.3) but clearer and with fewer excep- 
tions. 
2.1.3 The whole skull 

Twenty-seven variables were selected from the 
three anatomic regions to cover the morphology of the 
whole skull (Table 5). These variables were selected 
from those making the greatest significant contributions 
in the first two axes, whilst considering the balance be- 
tween the various regions. All PC1 eigenvectors are sig- 
nificantly positive, and the other axes include both posi- 


tive and negative eigenvectors. 


Table 5 Percentage of the total variation. eigenvalues and eigenvectors in the first two principal components based on 


the analyses of the whole skull 


Female Male 

1 2 I 2 
Eigenvalue 14,72 3.46 12.84 3.58 
Percentage 54.5 11.7 47.6 13.2 
Cum. per 54.5 66.2 47.6 60.8 
Eigenyectors 
LBCB 0.765 0.334 0.834 0.196 
LIAW 0.592 0.677 0.683 0.480 
CONMI 0.919 -0.197 0.881 -0.256 
MAMI 0.818 -0.293 0.675 -0.229 
UIIMDL 0.362 0.761 0.623 0.599 
UCMDL 0.835 0.070 0.743 0.166 
UIIBELL 0.581 0.414 0.509 0.603 
UMIBLLM 0.845 -0.111 0.701 -0.328 
LIIMDL 0.245 0.839 0.492 0.731 
LIAMDL 0.563 0.575 0.006 0.310 
LCBLL 0.569 — 0.006 0.636 0.239 
LMIBLLM 0.807 0.003 


The variables reaching significant contribution level are in bold. 













Female Male 
I 2 1 出 
MUZL 0.908 - 0.077 0.368 0.162 
BIORBW 0.743 -0.184 0.709 -0.382 
INTORBW 0.709 -0.118 0.559 -0.299 
PALWID 0.846 - 0.133 0.871 0.009 
UBCB 0.910 0.100 0.827 0.308 
UIAW 0.828 0.428 0.723 0.474 
PIRH 0.793 -0.335 0.758 - 6.302 
CALYL 0.508 - 0,253 0.852 ~ 06.355 
' MIDPARW 0.848 - 0.247 0.648 — 0.492 
INFRMAL 0.388 - 0,183 0.878 -0.251 
POSTORB 0.628 0.012 0.700 -0.149 
OCCH 0.489 -0.154 0.404 -0.239 
FORMAGL 0.469 -0.162 0.510 -0.423 
FORMAGW 0.806 - 0.249 0.260 - 0.336 
GLENOL 0.553 -0.105 0.492 0. 104 
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Females The total variance accounted by the 
first two axes is 66.2% (Table 5) .In broad terms the 
first two axes place three of the species in the third 
cluster ( M .thibetana , M . sylvanus and M . fuscata ) at 
the positive extreme, and all species in the first cluster 
( M. fascicularis , M . radiata and M . sinica ) at negative 
extreme (Fig.4).This picture resembles the first axis 
separations of the females in most of the functional units 
and overall regions (although it is less clear in the 
face) . Overall , this pattern is more like the masticatory 
region than the others. The second axis seems mainly to 
just separate individual female of species with no partic- 
ular pattern evident.The three clusters are significantly 
separated on the two axes. The variables showing signifi- 
cant positive influences to PC2 are LBCB, LIAW, 
UI1MDL, UILBLL, LIIMDL, and ULAW. All of these 
are anterior dental dimensions. In contrast to these 
eigenvectors is the height of the piriform aperture 
(PIRH) which has a significant negative eigenvector 
(> 0.300). 

@Males The first two axes accounted for 60.8% 
of the total variation (Table 5) .In broad terms the first 
axis contrasts with the picture in the females.There are 
almost no separations of clusters that resemble those in 
the various regions in the first axis (Fig.4).This pat- 
tern does not resemble particularly that from any single 
region (certainly, it is not,as is the case for females, 
similar to the masticatory apparatus ) . 

The second axis seems mainly to separate the 
species with no particular pattern evident. More vari- 
ables than in the females make significant contributions 
to this axis: Incisor alveolar breadths ( LIAW ), 
UI1MDL, UI1BLL, LIIMDL, LI2MDL, UBCB and 
UIAW are associated with the anterior teeth, and exhibit 
significant positive contributions. UM1BLL,LM1BLLM, 
BIORBW. PIRH, CALVL, midparietal breadth ( MID- 
PARW ) ,the length and breadth of the foramen magnum 
(FORMAGL and FORMAGW) have significant negative 
eigenvectors ; they are mainly associated with breadths of 
the cranium. 

@Both sexes 


sex according to overall size in a manner extremely sim- 


The first axis separates species by 


ilar to the masticatory region (Fig.5).The differences 


between the sexes lay in the second axis,a picture that 
is very similar to the corresponding axis in masticatory 
apparatus. This again replicates what is evident in the 
masticatory apparatus and is quite different from the 
other regions. However, it is in the remaining significant 
axes (third and fourth) that this analysis differs from 
that of the masticatory analysis, these remaining axes 
performing not inconsequential separations of many in- 
dividual species. 

3 Discussion 

Although three species clusters were first used 
heuristically for defining patterns in the analyses, they 
do show differences between functional units , anatomical 
regions and the whole skull. In addition, they all show 
major differences between the sexes. The variations in 
craniodental morphology between species or species 
group have frequently been considered to be associated 
with both functional adaptation and phylogenetic inerti- 
a,and resulting in taxonomic difference ( Albrecht, 
1976;Cheverud et af. , 1986; Harvey et al. ,1991;Hy- 
lander, 1979). 

Macaques are widely distributed from Africa to A- 
sia. They have adapted to a wide variety of ecological 
conditions and habitants. This high degree of adaptabili- 
ty can be related to a great behavioural and ecological 
plasticity . Macaques can thus serve as an ideal model to 
understand the relationship between habitat preferences 
and variables relating to their past and present evolu- 
tion, interspecific interaction ,and morphology . 

Craniodental variation, especially in masticatory 
apparatus between species might be relevant to adapta- 
tion to specific feeding bebaviour and dietary compo- 
nents .A significant difference in the masticatory appara- 
tus has been found referring to the primates with differ- 
ent dietary items, especially between the frugivorous and 
the folivorous species ( Hylander, 1979; Smith, 1983; 
Takahashi eż af ., 1994) .In order to determine whether 
the variation of cramodental structure between species 
clusters, especially that in the masticatory apparatus is 
mainly associated with the differences in functional 
adaptation or in phylogenetic inertia, it is necessary to 


examine the information available for the three species 
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groups. 
3.1 The first cluster of species 

The species in this cluster include M . fascicularis , 
M . radiata and M ,sinica. They are distributed in 
Southeastern Asia, India and Sri Lanka. Macaca radiata 
occurs in the low land and high land forests (Simonds, 
1965), scrub jungles semi-evergreen forest, bamboo for- 
est and rain forest ( Krishnan, 1971; Nolte, 1955) in 
south of India. Its populations frequently appear on the 
roadside near agricultural fields (Simonds, 1965 }. It 
feeds on wild fruits, flowers, leaves, seeds and other 
foods from agricultural fields and catch insects and 
small animals (Rahaman et al ., 1969; Krishnan, 1971; 
Nolte, 1955). In general, its adaptation in the wild is 
very similar to that of M.mulatta (DelClue et al., 
1992). Macaca sinica is an endemic species to Sri Lan- 
ka. It appears in all major forest zones in Sri Lanka, in- 
cluding arid, semi-evergreen transitional, and evergreen 
forests from near sea shore to 2 100 meters above the 
sea level (Fooden, 1979). The diet of this species is 
variable; they can process foods from as many as 46 
plent species (DeClue et al., 1992). Macaca fascicu- 
laris occurs in a wide variety of habitats in Southeast A- 
sia. Its distribution areas include the forests of tropical 
evergreen rain forests, riverine, mangrove, monsoon, 
bamboo , agriculture regions and other areas of scrub and 
grasslands (Eudey, 1980: Kurland, 1973, Southwick et 
al., 1972). This species utilises between 20 and 179 
plant species. Fruits are the major proportion (about 
51% —%5%),so that this species is a good seed dis- 
penser (Lucas, 1996). Leaves, stems, and bracts ac- 
count for a smaller food proportion (3% - 27%). In 
some cases insects and animal prey may account for as 
little as 4% and as much as 23% (DeClue et al., 
1992), 
3.2 The second cluster of species 

The species including in this cluster are M . mula- 
ta, M assamensis, M.nemesirina and M. nigra 
(when present in the data). These species occupy a 
large area in Asia. Macaca mulatta lives in a wide 
range of habitats than any other macaques. Its popula- 
tions occur both in forests with little human intervention 


and in highly urbanised environments. It is adapted to 
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deep snowfall, deserts, subtropical, tropical, sparsely 
vegetated habitats, both coniferous and deciduous 
forests, and swamps (Fooden, 1982; Southwick er al., 
1966,1968; Zhang et af.,1989).It is highly omnivo- 
rous (Lindburg, 1971; Makwana , 1979) . In some areas 
foods include 42.8% fruits,41.3% leaves and 15.9% 
flowers (Makwana, 1979), or 63% -70% fruits (Lind- 
burg, 1971 ). The distribution of M . assamensis includes 
the evergreen or dry evergreen forests, and bamboo, 
tropical and subtropical monsoon deciduous forests 
(Eudey, 1980; Fooden, 1971; Zhang et al ., 1989). The 
diet of M . assamensis contains mainly fruits, seeds, in- 
vertebrates and vertebrates. The main proportion in 
some areas, however, are leaves (DeClue et al. , 1992; 
Sarkar et al., 1996). The distribution of Macaca 
nemestrina roughly coincides with that of tropical 
broadleaf evergreen rain forests (Fooden, 1975) and 
secondary forests ( Southwick et al., 1972). This 
species is highly frugivorous (Fooden, 1971; Medway, 
1970;Rodman, 1978). At Menttoko , it eats fruits from 
22 plant species which are overlapped with those eaten 
by M fascicularis (Rodman, 1991}. Thus there is a 
great similarity in diet between M.nemestrina and 
M . fascicularis in some areas. M.nigra was seen to 
raid in maize field and in the tree of the primary forest 
{Groves , 1980). 
3.3 The third cluster of species 

The species in this group are distributed in North- 
em Africa { M .sylvanus ) and Asia (covering Northem 
India, Bangladesh, Burma, China, Malaysia, Thailand, 
Cambodia, Loas, Vietnam and Japan). The vegetation 
inhabited by M.sylvanus is a variety of grasses and 
herbs with little or no tree or shrub (Ménard et al., 
1990), mixed-oak forest (Fa, 1984, 1986), and rocky 
mountains and forests with a closed canopy and few 
shrubs (Fa, 1986; Menard et al., 1996). Foods from 
cedar trees can account for up 70% of the diet during 
times of deep snow, or deciduous leaves and herbs in 
other times (Drunker, 1984; Fa, 1984). The habitats of 
M . arctoides include primary the forests of evergreen, 
secondary evergreen, hill evergreen, mixed deciduous 
leaves, and dry rocky outcrops (Bertrand, 1969; Food- 
a, 1971; Treesucon, 1988). The diet of M . arctoides 
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include fruits, leaves and seeds from a variety of wild 
plant species and small animals { Bertrand , 1969; Food- 
en, 1971; Fooden et al. , 1985). According to McCann 
(1933) it is predominantly frugivorous. The habitats of 
M . thibetana are very similar to those of M . arctoides . 
They include subtropical broadleaf evergreen forests, 
temperate deciduous, and subtemperate mixed conifer- 
ous and deciduous forest { Fooden, 1983; Fooden et 
al.,1985).The diet of M . thibetana is composed of 
fruits , seeds , bamboo shoots and small animals { Xiong, 
1984; Zhao et al . , 1988). Macaca fuscata is distributed 
only in those parts of Japan covered with the deciduous- 
broad-leaved and subalpine conifer forest as well as 
bushes { Iwamoto, 1978; Izawa et al., 1963; Wada et 
af ., 1980). This species, however,is highly opportunis- 
tic in selecting its diet, fruits and seeds, over different 
times of the year and at different locations { Nakagawa, 
1989). 

From the information provided above, even though 
the habitats of macaques are variable and there are 
some preferences in diets referring to different species, 
there are no great variations in diet between the individ- 
ual species or clusters species found in this study. Nev- 
ertheless, the proportion of different food items in 
macaques are influenced by many factors, for instance, 
the diversity and availability of food resources, the re- 
quired energy to obtain such resources, and physiologi- 
cal and behavioural states {DeClue et af. , 1992). 

This basic commonality of diet implies that cranio- 
dental variation , especially in the masticatory apparatus, 
found between species clusters in this study may be less 
relevant to dietary adaptation and more associated with 
other factors,such as body size (Pan, 1998) .It is clear 
that the species in the first cluster are smallest and 
those in the third cluster are the largest ones. The 
species in the second cluster are intermediate. These 
vanations in body size among macaque species have 
been reported by many other researchers (e.g. Hill, 
1974; Napier et al . , 1967). 

It is not impossible that it is simply body size that 
might explain the main separation between clusters . The 
findings in this study indicate that M.radiata and 
M .sinica, of the sinica group, ate widely separated 


from M .arctoides , M . thibetana and M . assamensis in 
the same species group proposed by Delson { 1980). 
The former species are of smaller body size than the lat- 
ter due to divergent evolutionary pathways and adapta- 
tion to different geographical areas and environmental 
zones { Delson, 1980; Fa, 1989). Macaca radiata and 
M . sinica may have undergone dwarfing and as a result 
share few of their size-based characteristics with their 
putative group mates. 

However, the findings here in relation to differ- 
ences between the sexes {which are also largely related 
to differences in overall size in a data set that consists of 
measurements) imply that size is no a simple matter. 
For example , these two sets of size differences are large- 
ly orthogonal. Either one may appear in the first axis 
where all contributing variables are positive {this has 
been suggested to imply size), Both may appear (though 
oppositely arranged) in each of the first two axes. Fur- 
ther investigation of this size complexity is required. 

The relationships between some species revealed in 
this study are, however, similar to those from other re- 
search fields, The two stump-tailed macaque species, 
M . arctoides and M . thibetana ,are always closely clus- 
tered. The close relationship between them, as well as 
between them and M . assamensis {another species in 
the second cluster) ,has already been reported in stud- 
ies based on external morphology (Delson, 1980) , ge- 
netics (Zhang et ol . ,1990;Zhang, 1991 ) , biomolecular 
data (Cronin et al., 1980), craniometry ( Albrecht, 
1978), the similar incidence of lacrimal lacrimal fossa 
(Mouri, 1994) and hair structure { Inagaki, 1996) . This 
species is similar to M .arctoides and M . thibetana in 
this study. More details about this is available in Pan et 
al .(1998). 

Three other species, M . fascicularis , M . radiata 
and M . sinica ,form one cluster in each of skull analy- 
ses and show a great separation from the other two.This 
may imply that there exist close relationships among 
them from evolutionary and phylogenetic points of view 
as proposed by other studies. For instance, external fea- 
tures, especially the body weight and tail length { Pan, 
1998), gene frequencies (Melnick et al.. 1985; Weiss 
et al., 1973), relatively short genetic distances 
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(Cronin, et al., 1980), and similar allele-frequencies 
{Fooden et al. ,1989). 

Genetic studies show that M.nemestrina has a 
close relationship with the Sulawesi macaques including 
M nigra (Darga et al., 1975; Cronin et al., 1980; 
Melnick et al.. 1985). This study also demonstrates 
that M.nemestrina has a close relationship with 
M nigra being in the same cluster. Species related to 
M . nemestrina have been considered to be the ancestor 
of Sulawesi macaques. A Sundland population may have 
entered Sulawesi, resulting eventually in the macaque 
radiation in Sulawesi ( Delson, 1980; Fa, 1989; Fooden, 
1969, 1975). A genetic study supports this hypothesis 
(Kamamoto, 1996). 
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WH: 为 了 研究 狂 猴 属 的 颅骨 差异 性 ， 从 而 探 
讨 种 问 在 形态 、 功 能 和 系统 分 化 方面 的 相互 联系 ， 
测定 了 11 个 猕猴 种 类 的 77 个 颅骨 变量 ， 用 于 主 成 
分 分 析 和 判别 分 析 。 应 用 巢 式 分 析 方 法 ,分 折 过 程 
包括 3 个 步骤 。 所 有 变量 根据 功能 和 部 位 的 不 同 首 
先 分 为 了 个 单位 : Fa, Faw. km, be 
A. FRA, MASA. 2 步 根 据 它 们 所 
揭示 的 相似 性 (具有 相同 的 种 间 及 种 内 差异 性 类 
型 ) 合并 为 3 个 解剖 区 域 : RA CPA. FA 
Be. EG). EP EmA A) 和 整个 面 
AUG (HAUS AAR). 第 3 步 从 3 个 解剖 区 域 第 
选 出 27 个 变量 代表 整个 颅骨 的 形态 结构 。 除 了 各 
找 不 同 的 功能 单位 ， 解 剖 区 域 及 总 的 颅骨 具有 不 同 
的 种 间 和 种 内 差异 类 型 外 ， 此 过 程 对 证 出 研究 意义 
不 大 的 变量 起 很 重要 的 作用 。 上 述 分 析 过 程 分 别 用 
于 对 肉 、 寻 性 和 两 性 的 研究 。 所 研究 的 11 个 猕猴 
种 类 间 形 成 了 3 聚 类 。 第 1 RAR BR (Macaca 
foscicularis), BRE FR (M. sinica) ASK MR (M. 
radiata); 第 2 Rt (M. mudatta ), ARR 
(M. assensis), FHJ (M. nemestrina ) 和 黑 猿 
(M. nigra); 第 3 RUME (M. sylvanus), H 
Afk (M. fuscata), WEI (M. arctoides) 和 藏 
HI (M、thibetana }。 分 别 从 两 性 差异 、 食 物 、 生 


态 、 分 类 和 系统 分 化 方面 进行 了 差异 性 讨论 ， 结果 
认为 猕猴 种 间 颅 骨 的 差异 性 主要 是 由 于 系统 分 化 不 
同 而 引起 个 体 差异 所 致 ， 即 种 间 和 种 内 存在 的 个 体 
差异 。 在 主 成 分 分 析 中 ， 这 些 差异 在 不 同 的 区 域 表 
现在 不 同 的 成 分 上 。 在 咀嚼 器 官 上 种 间 的 差异 在 第 
1 主 成 分 上 上， 种 内 的 差异 则 在 第 2 ERSE. m 
的 情况 则 刚好 相反 。 这 两 种 差异 在 面 颇 后 及 颅 腔 上 
则 被 第 1 和 第 2 主 成 分 所 平分 。 这 样 ， 种 间 的 差异 
在 咀嚼 器 官 上 大 于 种 内 的 差异 。 种 内 的 差异 在 面 颅 
上 则 大 于 种 间 的 差异 。 这 两 种 差异 在 面 颅 后 和 颅 腔 
上 则 几乎 大 小 相等 。 这 一 研究 结果 表明 ， 与 传统 的 
概念 不 同 ， 第 2 主 成 分 不 仅仅 表现 形态 、 形 状 的 差 
F. 而 如 同 第 1 主 成 分 一 样 ， 也 表现 形态 的 大 小 成 
分 。 此 研究 所 揭示 的 猕猴 种 间 关 系 部 分 与 Foden 
(1976, 1980) 和 Delson (1980) 相同 。 如 平 顶 猴 
SER. SER. AERA ARM KA. PBR. 
头巾 猴 和 戴 帼 猴 的 关系 则 不 同 ,并 已 得 到 有 关 分 子 
生物 学 的 支持 ， 此 3 种 可 能 来 自 同一 祖先 并 经 历 相 
同 的 扩散 过 程 。 此 研究 所 设计 的 莫 式 分 析 过 程 提 殿 
了 一 种 很 好 的 差异 性 研究 手段 。 最 终结 果 暗 示 在 形 
态 学 研究 中 仅仅 考虑 某 一 区 域 的 形态 结构 是 很 不 够 
的 ， 因 为 不 同 的 部 分 具有 不 同 的 种 间 及 种 内 差异 类 
型 。 这 在 化 石 研究 中 尤其 要 注意 。 
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